Magnetosonic solitons in a dusty plasma slab 
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Abstract 

The existence of magnetosonic solitons in dusty plasmas is investigated. The nonlinear magne- 
tohydrodynamic equations for a warm dusty magnetoplasma are thus derived. A solution of the 
nonlinear equations is presented. It is shown that, due to the presence of dust, static structures 
are allowed. This is in sharp contrast to the formation of the so called shocklets in usual magne- 
toplasmas. A comparatively small number of dust particles can thus drastically alter the behavior 
of the nonlinear structures in magnetized plasmas. 

PACS numbers: 52.35.Bj, 52.35. Tc, 94.30.-d, 94.30Tz 
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In previous studies (e.g. Refs. [1, 2, 3]) of nonlinear magnetohydrodynamic waves, it 
has been shown that magnetosonic waves can appear both as sohtary waves and as shock 
waves in different kinds of plasmas [4-6] . The formation of magnetosonic shocklets in space 
plasmas has been discussed in recent papers [7]. For perturbations traveling across the 
external magnetic field direction, it was furthermore shown that exact analytical solutions 
of the magnetohydrodynamic (MHD) equations can be found for a cold plasma [5]. The 
nonlinear solutions have also been obtained for a warm magnetoplasma [6], as well as for 
investigating matter waves in dilute gases [8]. Later more general [9], or alternative [10], 
solutions were also presented. 

In the present paper we are going to deduce a nonlinear model for fast magnetosonic 
solitons (FMS) in a plasma where a layer of charged dust particles [11] is also present. It 
will then turn out that the presence of dust particles can significantly change the results 
of previous papers [2, 5] where electron-ion plasmas without dust were considered. Dusty 
plasmas can appear in many contexts, for example in interplanetary space, in cometary 
tails and comae, in planetary rings, in the Earth's atmosphere, as well as in dc and rf 
discharges, in plasma processing reactors, in solid-fuel combustion products, and in fusion 
plasma devices [12], and dusty layers [13] exist in the Earth's mesosphere |14| . 

The dynamics of the nonlinear FMS waves in a dusty magnetoplasma is governed by a 
set of equations composed of inertialess electron momentum equation 

= — ngcE — Vpe — — x B, (1) 

c 

the ion continuity equation 
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and the ion momentum equation 



+ ■ Vi = UicE - Vpi + riiC^ X B, (3) 



where ne{ni) is the electron (ion) number density, e is the magnitude of the electron charge, 
E is the wave electric field, B is the sum of the ambient and wave magnetic fields, Pe{Pi) = 
neTf^ijiiTi) is the electron pressure, T^iTi) is the electron (ion) temperature, Ve(vj) is the 
electron (ion) fluid velocity, nii is the ion mass, and c is the speed of light in vacuum. 
Following standard theory, we consider the plasma as isothermal, i.e. Tg and Tj are constants. 
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Equations (l)-(3) are supplemented by means of Ampere's law 

Ane , 



V X B 



(4) 



and Faraday's law 



dB 

'dt 
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(5) 



together with the quasi-neutrality condition Ue = rii — ZdUd. Here nd{x) is the prescribed 
number density of the dust, while Zd is the number of electrons on each dust grain. As the 
dust particles are assumes to be very heavy, we can consider them as immobile. CoUisional 
effects are also neglected. The quasi-neutrally condition holds for a dense plasma in which 
the ion plasma frequency is much larger than the ion gyrofrequency. Equation (4) is valid 
for FMS waves whose phase speed is much smaller than the speed of light. 
Eliminating E from (1) and (3) we obtain 11| 
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On the other hand, from (1), (4) and (5) we have 
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We are interested in studying the nonlinear properties of one- dimensional FMS waves 
across the external magnetic field direction zi?o, where z is the unit vector along the z 
axis and Bq is the strength of the ambient magnetic field. Thus, we have V = id/dx, 
Vj = u{x)^ + 'w{x)'y and B = B{x)z, where x (y) is the unit vector along the x (y) axis in 
Cartesian coordinates. Thus, with these prerequisites, the MHD equations become 
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where we have dropped the index i. We note that the magnetic field is frozen-in according 
to 

B/Bo = njno, (9) 

for some constant no, and that w —>■ as Ud ^ 0. However, these conditions are modified 
by the inclusion of dust species, leaving w ^ in the generic case. Thus, the inclusion of 
dust breaks the frozen-in-field line symmetry between the electron and ion species, similar 
to the effect of including inertial terms in a non-dusty plasma. As will be seen below, this 
has profound effects on the nonlinear dynamics of magnetized dusty plasmas. 

Next, we look for stationary solutions, i.e. dt = 0. With this, the above equations can be 
integrated. The continuity equation yields the velocity component u according to 

u{x) = ^, (10) 

n[x) 

where mq is some constant of integration determined by our boundary conditions. For the 
velocity component w we have the equation 

dw e ZdUdB 



dx mc {n — ZdUd) ' 
while from the momentum equation in the x-direction we obtain 
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Thus, the magnetic field is given by the quasi-neutrality condition and (9), the x-component 
of the velocity is given by (10), and it remains to solve (11) and (12) for n and w. We note 
that ffTTl) can be written in energy form according to dximw"^ /2) = {e/ c) ZdUdBw / {n — ZdUd) , 
and adding this to Eq. ([12]) gives dx[m{u'^ + w'^)/2 + T^lnne + Tlnn] + n-'^d^^B'^ /8n) = 0. 

Next we normalize our variables according to Ud ZdUd/nQ, n n/riQ, u —>■ u/uq, 
w ^ w/uo, Vt = {T/my/^ ^ Vt/uo, Va = {Bi/ATrmnoY/^ ^ Va/uo, ^ TJT, and 
X xuJc/uq, where Uc = eBQ/mc is the ion cyclotron frequency. 

Inserting (9) and (10) into (11) and integrating gives 



w{x) = / nd{x') dx' , (13) 



while integration of (12) gives 



1 
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FIG. 1: The normalized ion density n and magnetic field B plotted as a function of the normalized 
variable a; for a = 0.2, A^^ = 1, T/r = 0.1, Te = 10, and Va = 10. 
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FIG. 2: The normalized ion density n and magnetic field B plotted as a function of the normalized 
variable 2; for a = 0.2, Nd = I, Vr = 1, = 100, and Va = 10. 

where Wq is a constant. For a dust layer, we may use the density profile rid = 
(A''rf/v/27r a) exp(— x^/2a^), where A^^^ is the normalized number of dust particles and a is 
the width of the dust slab. Such a dusty plasma slab admits solitary wave solutions. In 
the figures, we show the numerical solutions to the normalized equations, using a gaussian 
dust density distribution. In Fig. 1, we have used a = 0.2, A^^^ = 1, Vr = 0.1, Te = 10, and 
Va = 10, in Fig. 2, we have used a = 0.2, A^'^ = 1, Vt = 1, = 100, and Va = 10, and 
in Fig. 3, we have used a = 0.2, A^^ = 1, Vr = 10, Te = 1, and Va = 10. In all cases, the 
velocity components stay bounded. 

To summarize, we have presented an analytical theory for magnetosonic solitons in a 
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FIG. 3: The normalized ion density n and magnetic field B plotted as a function of the normalized 
variable a; for a = 0.2, AT^ = 1, = 10, Te = 1, and Va = 10. 

magnetized electron-ion-dust plasma. It is found that the presence of a stationary charged 
dust layer provides the possibility of new classes of magnetosonic solitons. For typical 
dusty plasma parameters, we have displayed the density and magnetic field profiles of com- 
pressional magnetosonic solitons. The present results may describe the salient features of 
localized magnetosonic solitons in forthcoming laboratory experiments in strong magnetic 
fields. 
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